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New advances in experiments on the random-field Ising model, as realized in dilute antiferromagnets, 
have brought us much closer to a full characterization of the static and dynamic critical behavior 
of the unusual phase transition in three dimensions (d — 3). The most important experiments 
that have laid the ground work for our present understanding are reviewed. Comparisons of the 
data with Monte Carlo simulations of the d = 3 critical behavior are made. We review the current 
experimental understanding of the destroyed d = 2 transition and the experiments exploring the 
d = 2 metastability at low T. Connections to theories most relevant to the interpretations of all 
the experiments are discussed. 



1 Introduction 

The random-field Ising model (RFIM) has been an important focus of theoretical and 
experimental studies of the statistical physics of random and frustrated systems. Although 
there are some similarities, particularly at large random fields, to the physics of spin- 
glasses also covered in this book, the three dimensional (d = 3) ground state of the 
RFIM in the small random-field limit has the same long-range order as would be observed 
in the absence of random fields. Hence, the two models differ fundamentally. Nevertheless, 
the d = 3 RFIM transition is profoundly altered by the random field. For d = 2 the random 
field destroys the transition which takes place in the absence of the random field. Not only 
does the RFIM have significance in the formation of long-range order in real materials, 
where defects causing random fields are often present, it also challenges the methods and 
ideas of theorists and experimentalists that have been developed in past studies of phase 
transitions in pure, translationally invariant materials. There are a number of relevant 
reviews that have been written covering the formidable problems encountered in the exper- 
imental study of RFIM systems era. This one represents a comprehensive overview of the 
experimental situation in the most studied systems, the dilute anisotropic antiferromagnets, 
emphasizing the most current experimental results. The theories and computer simulations 
most relevant to interpretations of the behavior observed in dilute antiferromagnets will be 
included. A few systems that are not antiferromagnets will be mentioned in section 11. A 
comprehensive review of the theory ofthe RFIM by Nattermann also appears in this book. 

For d = 3 it has been rigorously B shown that a transition must take place for small 
random fields. As we shall see, the RFIM transition is very different from the more usual 
phase transitions encountered in antiferromagnets. The RFIM can be most simply modeled 
by spins on a lattice that point along one axis and are subjected to a random ordering 
field that competes with the long-range collective spin ordering. One simple Hamiltonian 
representing an Ising ferromagnet with an imposed random field is 

H = — JijSiSj — hiSi . (1) 

<i,j> i 
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The random field has the properties [h%] av = and \hf\ av = h% where [...] av denotes an 
average over the disorder. Most of the theoretical and simulation efforts, though not all, 
have focused on such ferromagnetic models. On the other hand, the most studied and best 
characterized experimental realization of the RFIM, by far, is the dilute, anisotropic antifer- 
romagnet in a uniform field applied along the spin ordering axis, which can be represented 
by the Hamiltonian 

H = X! Jij e i e jSiSj -^^HeiSi , (2) 

<i,j> i 

where ej = 1 if site i is occupied and if empty, and H is the uniform field. Locally, the 
sublattice with the most spins tends to align with the applied field in competition with 
long-range antiferromagnetic order in which one sublattice globally aligns with the field. 
The applied uniform field and the effective random field generated by it are proportionally. 
The random field is therefore easily controlled or even turned off completely. This provides 
the opportunity to do scaling studies not easily done in other systems. Importantly, samples 
can be cooled in zero field before applying the random field (ZFC). Other systems, such 
as those with structural phase transitions can only be cooled in the random field (FC). 
Since, as we shall see, hysteresis plays an important role in the understanding of the RFIM 
transition, the ZFC process is crucial. Of course, by virtue of critical behavior universality, 
the systems studied need not correspond precisely to the Hamiltonians above but must 
simply have the appropriate symmetries. 

Fishman and AharonyEl first noted that the dilute antiferromagnet in a uniform field is 
a RFIM system and Cardy§ showed that the critical behavior in the limit of small fields 
belongs to the same universality class as the uniform ferromagnet with random fields. These 
works opened up a tremendous opportunity to investigate the RFIM experimentally. An 
understanding of the RFIM phenomena in the dilute antiferromagnet is steadily evolving 
with experiments performed on very high quality anisotropic crystals. A major aim of this 
review is to present an overview of the d = 3 RFIM transition that takes place in dilute an- 
tiferromagnetic systems which is consistent with all of the published data (though certainly 
not all the published interpretations of the data). The d = 3 phase diagram has proven 
much richer than anticipated and this review necessarily encompasses high, intermediate 
and low magnetic concentrations as well as large and small random fields. The most recent 
experiments by Slanic, et al. H at high magnetic concentrations are promising as they ap- 
pear to afford the opportunity to make real headway in the experimental characterization of 
the RFIM critical behavior and in making connections to recent theoretical and simulation 
results. Such work is still in progress, so only preliminary results can be discussed. 

Theory and experiments on the RFIM have been closely tied throughout the period of 
investigation from the Fishman and AharonyQ work until the present, though there has not 
always been agreement. The greatest progress in the experimental investigations has come 
when a variety of techniques are employed and interpretations consistent with them are 
made. Often mistakes have been made when only one technique is relied upon for interpre- 
tation. A complication of studies using the dilute antiferromagnet is that random magnetic 
vacancies constitute strong pinning sites for domain walls a. Such strong vacancy pinning, 
while enriching the d = 3 antiferromagnetic phase diagram, is not present in the theoretically 
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well-studied ferromagnetic model. Random-field pinning, present in both antiferromagnets 
and ferromagnets, seems to be much weaker. The correspondence between antiferromagnets 
and the ferromagnetic models is best when the magnetic dilution is small, in which case 
the antiferromagnetic order is stable up to the transition, T C {H). For concentrations near 
x = 0.5, there is evidence that the long-range order breaks into static structure consisting 
of large, intertwined and weakly interacting domains well below T C (H). This has prevented, 
at these concentrations, a characterization of M s 2 vs. T and the line shape below T c (i£L 
When the percolation threshold is approached (x x v = 0.25) a de Almeida-Thouless Ii3 
behavior appears for larger H and the system appears to behave similarly to a spin-glass. 
In this review each of these three concentration regions is discussed. 

For the case of d = 2 dilute antiferromagnets, the random-exchange Ising model (REIM) 
transition is expected to be destroyedJj'O as soon as H, which generates the random field, is 
applied, and this has been observed EJ'EJ. The temperature regime well below the rounded 
transition, however, is still being investigated theoretically and experimentally 0. Both 
temperature regimes are briefly reviewed. 



2 Sample Considerations 

The most studied dilute d = 3 antiferromagnet suitable for RFIM studies is Fe x Zn\- x Fi. 
Its large crystal-field anisotropy persists Ii3 as the magnetic spins are diluted and it is there- 
fore an excellent Ising system for all ranges of magnetic concentration x. Crystals can be 
grown for all x with extremely small concentration variations 5x < 10~ 3 and with superb 
structural quality. The magnetic interactions are dominated by the second-nearest- neigh- 
bor super-exchange between the body-center and body-corner ions. All other interactions 
are negligible, except possibly near the percolation threshold concentration, where even tiny 
frustrating interactions become important£3li3. Another class of materials representing the 
anisotropic random-field systems is Fe x Mg\_ x Cl2- This system differs from Fe x Zn\- x F2 
in that it is layered. The layers are ferromagnetic and the interplanar antiferromagnetic 
coupling is comparable in strength to the intraplane coupling, making this a good d = 3 
Ising system. The smaller exchange in this system allows the large field region of the 
phase diagram to be exploredli^l. For x < 0.55, a strong second-nearest-neighbor competing 
exchange in the Fe x Mg\— x Cl2 system induces spin-glass behavior E3 and so random-field 
studies are restricted to higher x. There is excellent agreement between the random-field 
behavior of Fe x Mg\_ x Cl2 and that of Fe x Zn\- x F2- Some studies have also made use of 
the highly anisotropic Co x Zn\- x F2 system. A number of studies have been made in the 
less anisotropic system Mn x Zu\- x F2- The anisotropy in Mn x Zni- x F2, which is small for 
x = 1.0, decreases further upon dilution. Nevertheless, the H = REIM critical behavior of 
Mn x Zn\- x F2^ is quite consistent with that of Fe x Zn\— x F2^ and all of the RFIM experi- 
ments done on Mn x Zn\- x F2 seem qualitatively consistent with those done in Fe x Zn\_ x F2 
and Fe x Mg\_ x Cl2- The system does allow large applied fields relative to the anisotropy, al- 
lowing studies of the spin-flop regionclHa For the d = 2 RFIM, Rb2Co x Mg\_ x F± has been 
studied and appears to be an ideal system c2. It is very anisotropic and consists of layers 
of magnetic ions with a single dominant intralayer exchange interaction and an interlayer 
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interaction which is smaller by several orders of magnitude. 

Disagreements among the various interpretations of experimental data have arisen when 
concentration gradients obscured the true random-field behavior of a sample and were not 
fully appreciated in the data analyses. Although the gradient effects have been exten- 
sively reviewed c3'E3, the problem is still relevant to interpretations of recent experiments, 
as discussed below. Basically, one must realize that a concentration gradient will round 
a transition and can affect critical behavior measurements drastically. It is best if the 
gradients are unambiguously determined independently of the critical behavior measuring 
techniques. The size of the gradient can then be incorporated into the interpretation of 
the critical behavior data. Disagreements over interpretations of data in RFIM systems are 
usually resolved once the effect of concentration gradients are properly taken into account. 



3 Scaling Behavior Theory 

Although the scaling behavior of the RFIM has been discussed extensively in previous re- 
views, we emphasize the salient points again since many experiments are addressing the 
RFIM critical behavior and, unfortunately, not all of the current experimental interpreta- 
tions being proposed are consistent with scaling theory. Static critical behavior for temper- 
atures very close to the second-order transition temperature T c can generally be described 
by power law behaviors which become exact as the reduced temperature t = T/T c — 1 — > 0. 
We briefly list the ones most useful to us. The free energy has the asymptotic behavior 
|£| 2-Q , and the specific heat is correspondingly given by 



C p = A ± \t\~ a + B , (3) 

where we include a constant background term which describes the peak height when a < 0. 
For the case where a — > and A + j A~ — > 1, we use the symmetric logarithmic form 

C p = Aln\t\+B (4) 

instead. Several critical parameters can be obtained from neutron scatteringEl. The corre- 
lation length for antiferromagnetic fluctuations has the form 

-1 _ £±1 + 1-1/ 



S = K- L =g\t\-" • (5) 



The staggered susceptibility is 



xs = xt\tr . (6) 

For random-field systems we have the disconnected staggered susceptibility 

xs d = xo d± \tr . (7) 



The staggered magnetization, the order parameter for antiferromagnets, is given by 

M s = M Q \tf , (8) 
4 



which is only nonzero for t < 0. In these expressions + and — are for t > and < 0, 
respectively. The exponents and the ratios for amplitudes above and below T c are universal 
parameters common to all systems sharing the same symmetries. The asymptotic critical 
exponents satisfy scaling relations such as 

7 + a + 2(3 = 2 . (9) 

There are also hyperscaling relations that involve the dimension d such as 

a + vd = 2 (10) 

that hold for pure and REIM systems but are violated in the RFIM, in which case Eq. [l0| 
is modified by the additional violation-of-hyperscaling exponent, 0, with 

a + v(d-0) = 2 . (11) 

As the RFIM transition is approached by varying H or T, one observes a crossover from 
the zero-field universality class to the RFIM one. The crossover behavior can be described 
by a crossover scaling function. For example, the free energy is given by 

F = \t\ 2 - a f(th; 2 ^) , (12) 

where <f) is the crossover exponent, a is the zero-field exponent and h r is the random- 
field strength. A consequence of crossover between different universality class behaviors 
is that measurements may not yield asymptotic universal parameters unless the data are 
sufficiently close to T c . Rather, one obtains effective exponents. The scaling relations (not 
the hyperscaling ones) between exponents are still approximately satisfied by the effective 
exponents Ell. Another consequence of the crossover function is a depression of the phase 
transition temperature given by 

T C (H) = T N - AH 2 ^ - bH 2 , (13) 



where b represents a small mean-field shift also present in the pure system. The H-T 
phase boundary curvature is determined by 0. Fishman and Aharonyl showed that for the 
crossover from pure to random-field d = 3 behavior, (j) = 7, with 7== 1.25 obtained from 
theory and experiment H. Although some of the early experiments H'SEI were incorrectly 
interpreted as showing this, it was also arguedE3 that </> is much larger. The latter result now 
appears to be universal, with cf> = 1.42 + 0.03 obtained for Fe x Zni- x F<^^, (j) = 1.43 ± 0.03 
for Mn x Zn\- x F % ¥A and = 1.41 ± 0.05 for Fe x Mgi- x Cl 2 ^. Aharony @ predicted 
that for a random- exchange to random-field crossover, eft is several percent larger than 7. 
Neutron scattering measurementsE3 in -Feo. 46 ^^0.54-^2 yielding 7 = 1.31 ±0.03 confirm this. 
This was similarly verified Ej in Mn x Zn\- x F2 with the result 7 = 1.36 ± 0.08. The early 
interpretations EUSEl that <fi = 1.25 were influenced by the concentration gradients in the 
samples c3 and the resulting misidentification of T/y. 

The scaling function has other consequences for random- field antiferromagnets. Fish- 
man and Aharony Q obtained 

kT x = A 1 +A 2 \t\ l ~ a -A?,\t\ 2 P (14) 
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for the static H = uniform susceptibility for a system dominated by the pure critical 
exponents /3, a and 7 at H = 0. The experimental systems, however, are dominated by 
random- exchange critical behavior at H = and the correct relationship is therefore @ 

kT X = + A 2 \t\ 1 ~ a - A 3 \t\ 2 - a ~' t> , (15) 

using the REIM exponent a and the REIM to RFIM crossover exponent 4>. For H > 0, 
Kleemann, et al. showed, by considering leading singularities in derivatives of the free 
energy for H 7^ 0, that the field dependence of the amplitude of the peaks in the specific 
heat, (dM / 8T)h and (dM/dH)x is governed by the exponents y = (2/</>)(a — a), y = 
(2/(/>)(l + a — a — 4>/2), and y = (2/c/>)(2 + a — a — </>), respectively, where a is the RFIM 
specific hea t exponent. The exponents can be determined from specific heati and Faraday 
rotation e3 experiments on Fe x Zn\- x F2 to be y « 0.1, y = 0.60 ± 0.10, and y rs 0.97, 
respectively. For Feo.7Mgo.3Ch, y = 0.56 ± 0.05§ for (8M/dH) T . Since the exponents 
(j) = 1.42 ± 0.03 and a = —0.09 ± 0.03 are known fairly accurately, we can invoke these 
results as a strong indication from scaling that a ~ in good agreement with direct specific 
heat results discussed below. There is ample evidence that scaling works well in all of these 
systems, despite recent arguments £3 to the contrary. Hence, at this time, experimental 
interpretations should be constrained to agree with scaling theory. 



4 The d = 3 RFIM Transition 



The first evidence that a phase transition occurs in the d = 3 Ising model came from the 
critical behavior of the specific heat measured using optical linear birefringence H (An) 
techniques, which minimize the effects of concentration gradients since the laser beam is 
directed perpendicularly to the concentration gradient. The proportionality ElJll between 
and the magnetic contribution (C m ) to the pulsed specific heat (C p ) data has been 
shown explicitly ljQ for Feo.4eZno.54F2 and Feo.93Zno.07F2- In anisotropic systems the 
Zeeman contribution to the birefringence is negligible at reasonable fields e3. In addition, 
Faraday rotationH and susceptibilityS measurements yield the sp ecific heat critical behav- 
ior. Recent claims to the contrary eH'e3 have no theoretical basisc3 and depend on analyses 
of experimental data which have been questioned e^I. 

The specific heat critical behavior in the intermediate range, 0.4 < x < 0.8, has been 
measured in Fe x Zn±- x F2 @c3 and Fe x Mg\_ x Cl2 c3 with enough precision to determine 
that ZFC data are well described by a symmetric logarithmic divergence over a reasonable 
range in t. At very small \t\, rounding is observed and is attributed to the tremendous 
critical slowing down of the RFIM in dilute antiferromagnets, as will be discussed in the 
following section. FC always yields behavior that is much more rounded because finite- 
size metastable clusters @S are frozen in just above T C (H). The region over which the 
logarithmic ZFC behavior and the dynamic rounding are observed increases with H as 
expected from crossover scaling (Eq. 12) and dynamic scaling as discussed in the next 
section. 

Figure 1 shows recent and C m data for a high magnetic concentration sample a 

Feo.93Zno.07F2- We have subtracted the approximate large phonon contribution to C p , leav- 
ing only C m , so that the correspondence of the two techniques is apparent. We accomplished 
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Figure 1: d ^j£ > vs. T, where A — 9.17 x 10 -6 K _1 is the same proportionality constant found for pure FeF2, 
and C m / R vs. T for ^0.93^710.07^2- The specific heat has the phonon contribution subtracted as discussed 
in the text. ZFC data are shown in the main figures. The insets show the H = 7 T FC data as well as the 
ZFC data for the case. The curves in the left figure are the same as the curves in the figure on the 

right except that they are rounded by the larger, measured concentration gradient. For the specific heat 
inset, the ZFC data are not shown, for clarity, but the dotted line is the same as the solid ZFC line in the 
main figure. Just as in experiments at lower concentrations, the birefringence and pulsed heat techniques 
yield precisely the same behavior, including FC-ZFC hysteresis very close to T C (H). The critical behavior 
for H > is closely approximated by a symmetric, logarithmic divergence. 



this by assuming the same proportionality between C m and ^ found for both Fei 7 ^ 
and i ? eo.46^o.54^ 7 2- The excess specific heat contribution found for the H = was then 
subtracted from all the C p data sets. The solid curves in the C m figure are adapted from 
those in the -^y-^ one by first numerically rounding the data by the known gradient, 

drawing smooth curves through them, and then transferring the curves to the C m figure 
with no other adjustments. The C m peaks are not as sharp since the entire sample was 
used, increasing the effects of concentration gradients relative to the data which are 

sensitive only to the gradient along the laser beam. Clearly, the known gradient accounts 
well for the difference in the C m and jy data. The insets for both sets of data show the 
FC behavior at H = 7 T. The FC data are more rounded than the ZFC. The FC C m data 
are shown in the inset on the right with the FC curve, rounded by the concentration 

gradient in the same manner as the ZFC ones, shown as a solid curve. The FC curve corre- 
sponds well with the C m data. The dotted curve corresponds to the ZFC and is the same as 
the solid one in the main figure. For clarity, we do not show the ZFC data in the inset. No 
other adjustments have been made. Clearly, the hysteresis is much more difficult to discern 
in the specific heat datac3, but this is consistent with the larger concentration gradient. The 
hysteresis can only be observed in samples with extremely small gradients. This is certainly 
one reason why some experiments on samples with appreciable concentration gradients fail 
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t the and C m data yield the same 

for Feo.46 Zno.5iF2. As discussed below, 



to exhibit hysteresis. In all respects we see t 
critical behavior just as was found previously 
-Feo. 93^^0.07-^2 yields neutron scattering line shapes that are fundamentally different from 
those obtained at lower magnetic concentrations in that they do not show hysteresis well 
below T C (H). Yet, the specific heat appears remarkably similar to that of lower magnetic 
concentrations. The symmetric, logarithmic (a = 0) behavior for ZFC contrastswith the 
H = data that are consistent with the expected asymmetric random-exchange Eil cusp. 

For all concentrations there is a temperature, T eq (H), below which hysteresis between 
the FC and ZFC procedures plays a role in the specific heat as well as most other ex- 
periments. We emphasize that C m and data show precisely the same hysteresis, 
contradicting recent claims cH. Using a capacitance technique E3 on Feo.46 -2^0.54-^2 and 
Fcq, 72^^0.28-^2) the equilibrium boundary T eq (H) has been shown to lie just above T C (H), 
scaling precisely in the same manner with eft ~ 1.42. The nature of this boundary is still not 
entirely clear, though it is sharp enough to be measured precisely. It could be related to the 
extreme critical dynamics discussed in the next section or it could be related to a RFIM 
spin-glass-like behavior between T eq {H) and T C (H). Certainly one must be careful 
about the data extremely close to T C {H) since the system could be out of equilibrium. 

Hysteresis in the specific heat is not well observed in ac techniques il used on the 
less anisotropic system Mn x Zn\- x F2 at 6.6 Hz. The extremely rapid heating and cooling 
method (1QK per minute) of measuring specific heat in Fe x Mg\_ x Cl2 also shows very little 
hysteresis S, although early neutron scattering!^! and measurements in this system clearly 
exhibit hysteresis. Perhaps the time dependent techniques obscure the difference between 
FC and ZFC, though this is not yet clear. Recently, it was claimed that in Feo. 5^/10.5-^2 n o 
hysteresis is observed in the specific heat0. Although there is no published description of the 
procedures used, some conjectures can be made as to why the hysteresis was missed. Perhaps 
the sample concentration gradient induced rounding of 0.3 Kc3 obscures the transition 
the low field c3. The phase boundary might have been exceeded at the high field 
Finally, if the measurements were not sufficiently adiabatic, th e hys teresis may be obscured 
as they appear to be in other time dependent measurements l 55 H 34 l. The answer is simply 
unclear at this time and the failure to observe hysteresis could be a combination of effects. 

Whereas all the high resolution specific heat measurements done to date indicate a 
symmetric, logarithmic divergence with no evidence for any accompanying background dis- 
continuityil, Monte-Carlo simulationsil'iy indicate a cusp, with a large, negative exponent. 
This discrepancy between the exponents from simulation and experiment is, as yet, unre- 
solved and is certainly a major challenge to be addressed. 

In contrast to the birefringence measurements that first showed evidence B@ of a d=3 
transition, early neutron scattering measurements obtained with the FC process were in- 
terpreted as indicating a destroyed transition £3. Upon FC, no Bragg scattering is observed 
for concentrations x < 0.8. Instead, a finite-width shape approximated by a squared- 
Lorentzian, as discussed below, appears. We now know that long-range antiferromagnetic 
order is difficult to establish upon FC at low concentraton, but that a phase transition is 
nevertheless the basis of the underlying physics. Long-range order does occur for T < T C {H) 
when the field is applied after ZFC and FC domains have been shown to be metastableS. 
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Quite different phenomena are observed at high magnetic concentration. Recent scat- 
tering measurements cl using a crystal of Fe$. 93^770.07-^2 indicate a ZFC transition that is 
as sharp as allowed by the concentration gradient 5x = 0.002. More importantly, there is 
no evidence of nonequilibrium hysteresis except extremely close to T C (H), as in the specific 
heat Ex What is most remarkable is that the Feo. 93-^770.07-^2 neutron scattering line shapes 
show little hysteresis at low T. For x < 0.8, such hysteresis has always been observed and 
has been a major obstacle to interpreting the critical scattering below T C {H). An important 
distinction can be made between hysteresis seen in line shapes at lower concentration well 
below T C {H) which most likely originates in the large number of vacancies, and the rounding 
near T C {H) that appears even at high magnetic concentration. The latter may well have to 
do with RFIM critical dynamics as is clearly the case with the specific heat behavior cB. 

The interpretation of the scattering results in RFIM studies is severely hampered by 
the lack of adequately characterized line shapes provided by theory. As previously reviewed 
in more detail Bo, mean-field theory yields an elastic scattering cross section of the form 



S(q) = X(q) + M s 2 5(q) = 3— ^ + M s 2 5(q), (16) 



A 
q z + k 

for a pure system and, with an additional squared-Lorentzian term, 



S(q) = X(q) + X d (q) = ^2 + {q 2 + K 2)2 + M Ml), ( 17 ) 

for a random-field system HI. These expressions can be only approximate for d = 2 or 
d = 3 in pure or random systems, as one can see from the required asymptotic behaviors 
x(0) ~ K 2_r ' ~ \t\~ 7 and x d (fy ~ M s 2 ~ ~ \t\~^- The correspondence between the 

measured line shapes and the mean-field line shapes in pure systems is fairly good for 
d = 3 since 77 ~ 0.Q4 is small, though evidence for deviations from mean-field behavior 
have been observed e3. Although Pelcovits and AharonyE^I predict significant deviations 
from the Lorentzian line shape for T < T C {H) in the d = 3 REIM, where rj is also small, no 
definitive evidence for this has yet been observed in experiments. For d = 2, the discrepancy 
between the line shapes of the pure system and mean-field theory is more evident l3'E3 since 
i] = 1/4. For random-field systems, 77 ~ 1/2 is large S and the mean- field terms in Eq. 17 are 
expected to be far from accurate. The observed line shape in the random-field systems is in 
many cases inconsistent with the Lorentzian in Eq. 16, as was first shown by Yoshizawa, et 
al.lM However, the story is not as simple as adopting Eq. [l?] since this expression is often 
inconsistent with the dataS, particularly below T C (H). Nevertheless, Eq. 17 is a start. 



The first neutron scattering observations of the d = 3 RFIM phase transitioned, made 
using Feo.6-^^0.4-^2, showed that above the transition the ZFC line shapes are inconsistent 
with the single Lorentzian term but can be fit adequately by the sum of Lorentzian and 



squared-Lorentzian terms as in Eq. 17. Non-Lorentzian line shapes had been observed 
previously ti3 for the metastable domain state well below T C {H) after FC. From the fits to 
the critical behavior above T eq (H), the estimations v = 1.00 ±0.03, 7 = 1.75 ±0.20 and 7 = 
3.5 ±0.3 are obtained^. An attempt at a better characterization, including for T < T C {H), 
was madeS with the very homogeneous crystal i ? eo.46-^rao.54-F2. Several scans are shown 
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Figure 2: Neutron scattering intensity, I(q) vs. q for Feo.iaZno.54F2 and F eo.93Zno.07F2 after ZFC. Above 
T C (H), determined from the peak in the critical scattering, both samples exhibit Lorentzian plus squared- 
Lorentzian line shapes. Below T C (H), the Feo.46Zno.54F2 crystal shows a resolution limited lineshape which 
conforms to neither a Lorentzian nor a squared-Lorentzian line shape in addition to the Bragg scattering 
peak. The Bragg scattering decreases dramatically at this concentration as T C (H) is approached and shows a 
large hysteresis upon temperature cycling below T C (H). This behavior is consistent with the system breaking 
into large, intertwined, weakly interacting domains, a result of the very large number of vacancies at this 
concentration. In contrast, the Feo.g3Zno.07F2 crystal exhibits Lorentzian line shapes below T C (H) with no 
sign of domain formation. No hysteresis for \q\ > is observed well below T C (H), indicating equilibrium 
behavior. The Bragg intensity remains large extremely close to T C (H), indicating that /3 is very small, in 

agreement with simulations. 



in Fig. 2. Although the ZFC scattering above T C {H) is indeed fairly well fit by Eq. 17, the 
scattering line shapes below T C {H) are certainly not. Below T C (H) the measured scattering 
profiles are much too narrow, being essentially resolution limited for all T < T C {H) instead 
of having a width that varies as k(T). Furthermore, the intensity of the Bragg component is 
surprisingly small near T C (H) and qualitatively it appears as if the Bragg peak transforms 
into the non-Lorentzian, resolution- limited scattering contribution as T C {H) is approached 
from below. Finally, a peak in the q = scattering intensity is observed H at T C (H) upon 
ZFC and the height of the peak grows approximately logarithmically with time, a result of 
the extremely slow dynamics. Such a peak is difficult to observe normally but is evident in 
this case because the Bragg scattering, which usually swamps the critical fluctuation peak, 
is abnormally small just below T C {H). 

Although it is clear that the scattering is peculiar and interesting, extinction effects 
give pause to direct interpretations of the Bragg scattering intensities in the ^60.46-^^0.54^2 
experiments. High quality crystals scatter neutrons that are precisely aligned for the Bragg 
scattering condition in the first ten microns or so of material. As the scattering cross 
section diminishes upon approaching T C (H), the scattering simply occurs over a larger 
volume. Hence, the scattered intensity is saturated and does not exhibit the power law 
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behavior in Eq. || This difficult problem has been overcome by examining!!!! an epitaxial O 
Fe$. 52^^0.48-^2 fihn of thickness 3.4 ^m, grown on a ZnFi substrate. The film is thin 
enough to avoid extinction effects but thick enoug h (w 10 4 lattice spacings) for d = 3 
critical behavior. 

The neutron scattering results for the H = Bragg intensity^ of the film are consistent 
with REIM behavior. Hence, the film is high enoug hEi in quality to reliably reflect the 
d = 3 critical behavior. The scattering intensity for q > 0, coming solely from a Lorentzian 
contribution, is too weak to be observed in the film. For H > the scattering results are 
highly unusual. The ZFC Bragg intensity vs. T has the opposite curvature to that observed 
for H = 0, so the Bragg scattering intensity is very small quite far below T C (H). The loss 
of the ZFC Bragg intensity is irreversible below T C (H). This behavior has been interpreted 
as the system breaking into two intertwined domains with equal numbers of spins in a 
similar pattern to that observed in FC simulations El at low T. The formation of domains 
is observed to be irreversible below T C (H), a result that is consistent with the irreversibility 
observed in magnetization and optical studie£3c3 E. Inside the domains the spins are well 
ordered. The domain walls at this concentration (x ~ 0.5) are able to pass predominantly 
through the numerous vacancies, costing the system very little energy. It is clear that the 
Imry-Ma domain wall energy arguments S fail here since the energy needed to create such 
a domain wall is insignificant compared to the Zeemann energy decrease. Furthermore, the 
domains are only weakly interacting and each contributes to the phase transition at T C (H). 
Since the domains form well below T C (H), neutron scattering measurements are unable at 
this concentration to determine the critical behavior of the order parameter. Another piece 
of evidence indicating that the hysteresis for T C {H) comes from domain formation is found 
in the experimental results of x-ray scattering studiesE! at the surface of Mno.75.Zno. 25-^2- 
In the presence of surface defects, no hysteresis is observed, most likely a result of the defects 
preventing the formation of the two intertwining domains. When an identical sample was 
polished, removing the majority of defects, the hysteresis reappeared. 

In both the film (Feo.52 -Zno.48-^2) and bulk (-Feo.46-Zno.54-^2) studies, we find a large 
resolution-limited scattering line shape below the transition that is not well fit by either a 
Lorentzian or squared-Lorentzian term. It is most likely that this non-Lorentzian scattering 
profile is a signature of domain structure that forms below T C {H) even upon heating after 
ZFC. With this structure present it is very difficult to determine the critical behavior of 
the RFIM below the transition. Local probes like NMRl9, Mossbauer or /iSR in principle 
could yield the order parameter critical behavior, but prove to be complicated because of 
the spatial variations within the system. This motivated an investigation at a much higher 
concentration, where the vacancy concentration is small enough that domain walls cannot 
easily avoid a large energy cost of formation. Since hysteresis at low temperatures is seen in 
the work0 on M no. 75-^^0.25-^2, it is clear that one must go to even higher concentrations. 
Preliminary measurements 19 using the -Feo. 93^^0.07-^2 crystal seem to confirm the idea; the 
hysteresis in the scattering profile at low temperatures is eliminated. 

The abrupt change in line shape of Feo.93^no.o7-^2 at T C (H) is striking. Figure 2 shows 
scans taken just above and just below T C (H) at H = 7 T. Just 0.13K below T C (H) the line 
shape is incompatible with any significant squared-Lorentzian term. A Lorentzian term fits 
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Figure 3: k vs. T and \ vs - T for Feo.93Zno.07F2 for H = and 7 T, obtained from preliminary analysis 
of the neutron scattering line profiles for \q\ > 0. For H = 7 T and T > T C (H), the |g| > scattering is fit 
to a Lorentzian plus squared-Lorentzian lineshape. A Lorentzian was used in all other cases. For H = 7 T, 
the open triangles are for ZFC and the filled ones are for FC. The lack of hysteresis indicates equilibrium 
behavior. The solid curves are fits to the data. However, for H = 7 T and T < T C (H), no power law 
describes the data well, so no curve is shown for this case. 



fairly well. The absence of the non-Lorentzian component is most likely a signature of the 
stability of the long-range order right up to T C (H). Above T C (H), on the other hand, the 
line shapes are much more compatible with a fit to a Lorentzian plus squared-Lorentzian as 
in Eq. 17. The abrupt disappearance of the Bragg peak at T C (H), indicating a very small 
value for 0, contrasts greatly the behavior observed at lower concentrations in Fe x Zri\- x Fi. 
A small value of is consistent with theory and simulation resultslHSEl. The only previous 
experimental measurement El of is from dilation experiments on the lower concentration 
sample Feo.46ZriQ.54F2 which indicates < 1/8. This suggests that the small exponent 
value holds for lower concentrations even though the neutron scattering Bragg intensity 
cannot show it. The small value of is perhaps suggestive of a first-order transition >_but 
no latent heat is observed in the specific heat in the experiments or simulations i@@E§. 

At the lower magnetic concentrations, severe hysteresis is observed in the line shapes 
below T C (H). In the case of Feo.Q3Zrio.07F2, however, the line shapes for q > do not 
exhibit hysteresis except for the region near T C (H) where critical dynamics dominate. The 
Bragg intensity does show some hysteresis, being somewhat larger upon FC, but this is an 
extinction effect reflecting the fact that long-range order on length scales well beyond the 
instrumental resolution is not established upon FC, most likely a result of RFIM dynamics 
very close to T C (H). The temperature dependence of the Bragg intensity is essentially the 
same for the Bragg intensity upon ZFC and FC well below T C (H). 

Evidently, if we can extract the critical behavior in the Feo.93Zrio.oYF2 sample, it should 
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represent the equilibrium behavior since it is history independent. Unfortunately, the lack 
of a theoretical scattering line shape that goes beyond the misleading mean-field theory of 
Eq. 17 has severely limited the extraction of critical parameters below T C (H). Experimental 
work in this area is ongoing with progress anticipated, but theoretical work is also much 
needed in the near future. Above the transition the fits to Eq. 17 seem to work fairly well 
and one can extract the exponents, albeit with trepidation regarding exact results. The 
results for k and x vs - T are shown in Fig. 3 along with fits represented by the solid 
curves. Fits were made for all of the H = data and for T > T C (H) with the H = 7 T 
data. No suitable fit to a power law is obtained for T < T C {H) and no curves are shown. 
Preliminary fits! for T > T C (H) yield v = 0.93 ± 0.03, 7 = 1.71 ± 0.06 and 7 = 3.0 ± 0.1 
for 10~ 3 < t < 10~ 2 . These values are in reasonable agreement with earlier experimental 
resultsU at x = 0.6 mentioned above but are in disagreement with other estimations where 



the transition appears distinctly rounded Ell from concentration gradients. (Larger values 
for v have been obtained in other studies, but only because T C {H) has been taken to be well 
below the minimum in k in samples with relatively large gradient induced rounding.) There 
is reasonably good agreement between the exponents obtained from neutron scattering 
-Feo. 93-^710.07-^2 and those obtained from Monte Carlo simulations. For example, Rieger 
obtains v = 1.1 ± 0.2, 7 = 1.7 ± 0.2, 7 = 3.3 ± 0.6, and = 0.00 ± 0.05 for a Gaussian 
distribution of random fields. Thescattering results are also reasonably consistent with 
recent high temperature expansion Ea results for 7 and 7. 

Keeping in mind the uncertainty concerning the scattering line shape appropriate for 
analyzing the FeQ.93Zno.07F2 data, the preliminary scattering exponents above T C (H), 
7 = 1.71 and and from the specific heat, a ~ 0, satisfy the simple scaling relation in 
Eq. H if (3 is small as expected from theory. In stark contrast, a typical result from Monte 
Carlo simulations is that a is large and negative, for example a = —0.5 db 0.2 £3. Nev- 
ertheless, the specific heat exponent is the most consistent experimental exponent. Note 
that the measured amplitude ratio A + /A~ is very close to unity which is consistent with a 
logarithmic divergence. Also, as demonstrated in section 3, amplitude scaling relations for 
dilute antiferromagnets strongly indicate a ~ 0. 

We can use Eq. 11 and the measured exponent v = 0.93 for T > T C (H) to estimate 
the violation of hyperscaling exponent 9 = 0.85. Using the relations 7 = v(2 — rj) and 
7 = v(A — fj) with the values from scattering 7 = 1.71 and 7 = 3.0, we can estimate 
rj = 0.16 and fj = 0.77. These values are smaller than theoretical estimates, but they are 
very preliminary and further measurements and analysis will certainly refine them in the 
near future. The point to be made is that we are finally almost at the stage where serious 
comparison with theory can be made, though we are greatly hampered by not knowing the 
correct line shape. 

Finally, we should briefly mention a very recent suggestion by Birgeneau, et al. 

EH 

that 

the unusual curvature of the Bragg intensity versus T is actually a rounding of the phase 
transition at intermediate concentrations - the "trompe l'oeil" phenomenological model, as 
they have labelled it. It was introduced in an attempt to describe the scattering, mag- 
netization and specific heat behavior of the d = 3 RFIM phase transition in the lower 
concentration antiferromagnets Feo.5Zno.5F2 and Mno.75Zno.25F2- The interpretation of 
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Figure 4: x' vs. T at / = 1 Hz for Feo.47Zno.53F2 for applied fields H = 0.8, 1.6, 2.4, 3.2 and 4 MA/m. 
Open circles are ZFC data and filled circles are FC. The main contribution to the peak is from Griffiths- like 
clusters which form above T C (H). The smaller peak, which is resolved only at larger fields and only upon 
FC, is at the actual phase transition. From the dependence of the small ZFC peak height vs. the frequency, 
the dynamics can be ascertained. The behavior is consistent with a power law with a very large dynamic 
exponent or with activated dynamics. The inset shows a fit to the broad peak at H = 3.2 MA/m using a 

phenomenological Griffiths-cluster model. 



the data in this model conflicts sharply with the interpretations presented in this review^, 
since it clearly violates scaling for H > 0, which was developed by Kleemann, et al.E^ and 
Fishman and Aharonyi and is described in section 3. The authors take this as evidence 
that the scaling theory is incorrect. The interpretation requires that the peak in {dM / dT)jj 
coincide with the peak in (dM s 2 / 'dT)u and, to accomodate this, the (dM s 2 / dT)n data are 
adjusted within the thermometry uncertainties. The shifts of the data weaken the motiva- 
tion for the new model and the argument that scaling fails. The proposed model requires 
that the uniform magnetization couple strongly tothe antiferromagnetic long-range order 
and this has not yet found theoretical motivation e3. The model also depends on the spe- 
cific heat in Feo.5Zno.5F2 showing no hysteresis^ but such hysteresis has been observed in 
pulsed heat experiments using Feo.4eZno.54F2 and Feo.93Zno.07F2 S crystals with very 
small concentration gradients. 



5 Critical Dynamics of the d = 3 RFIM Transition 



The critical dynamics of the d = 3 RFIM transition are extraordinarily slow. Many of 
the experiments in Fe x Zn\- x Fi for H > that would normally be considered static mea- 
surements have shown behavior with approximately logarithmic_time dependence. These 
include neutron critical scatteringtia and capacitance techniques^. Spin-echo neutron scat- 
tering techniques showEJ that very small fields suffice to freeze the system over the entire 
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critical region \t\ < 0.1 in the nanosecond time regime. The most direct measurements 
of the RFIM critical dynamics are of the peak height of the ac susceptibility. The first 
susceptibility measurements on a RFIM antiferromagnet were performed on GdAlO^ : La 
by Rohrerta. Although at the time it was thought that the very rounded transition was 
evidence of the destruction of the d = 3 RFIM transition, it is now understood that the 
rounding is caused by slow dynamics. King, et al. measured the peak height in the ac 
susceptibility of -Feo. 46 -^"-0.54-^2 as a function of frequency and showed that the behavior is 
consistent with either a power law behavior 

x'(u;)~\t\- a F(u;\tn , (18) 

where zv ~ 14 has an unusually large value, or with activated dynamics with 

XV) ~ \t\- a GQnJ) , (19) 

where 9 is the violation-of-h yp erscaling exponent (Eq. |TT| ), as predicted by Villain il and 
Fisher 0. Later Nash, et al.H extended the measurements on the same sample to a very 
large frequency range of 5 x 10~ 3 < lu/2tt < 10 5 Hz and showed that activated dynamics are 
favored by the data with 9 = 1.05 ± 0.2. This is in accord with the violation of hyperscaling 
relation {d — 6)v = 2 — a using the measured values of v ~ 1 and a ~ 0, though it has 
been suggested that corrections to scaling should be considered 0. However, the picture 
changed substantially when, recently, Binek, Kuttler and Kleemann@ demonstrated that 
in -Feo.47 -2^0.53-^2 the peak in the ac susceptibility studied previously is not that of the 
phase transition itself but rather is due primarily to the dynamics of Griffiths-like spatial 
fluctuations E3 above T C (H). It was shown that the true critical peak corresponding to the 
phase transition is but a small peak that is not resolved at low fields and was therefore 
missed in earlier studies, as shown in Fig. 4. The shape of the peak is consistent with 
the exponent a ~ obtained in other experiments. High resolution measurements for a 
frequency range 3 x 10 _1 < u < 3 x 10 3 Hz again show that the peak is adequately described 
by either the power law with an unusually large exponent, zv ~ 14, or with activated 
dynamics. The critical peak, visible only upon ZFC, is surprisingly weak, indicating that 
only a small portion of the spins are involved in the phase transition. This is consistent 
with the very small peak observed in specific heat experiments at this concentration @. 
The larger peak has been related c3 to Griffiths-like instabilities in the temperature range 
between T C {H) and T/v, as discussed in the next section. Further refinement of the theory for 
the ac susceptibility peaks and investigation of other samples, for example -Feo.93Zno.07F2, 
may eventually settle the question of which dynamic model best fits the d = 3 RFIM 
in dilute antiferromagnets. The unusual RFIM dynamics have also been observed e3 in 
-Feo.7Mgo.3Ci2 using Faraday rotation techniques, where a symmetric logarithmic peak is 
seen with rounding. A fit of the peak height to a power law behavior yields zv = 8.3. 

6 The d = 2 Destroyed RFIM Transition 

In contrast to the d = 3 case, it is clear from theory 00 that the d = 2 phase transi- 
tion is destroyed by the random field. Experimentally this was demonstrated definitively 
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Figure 5: d(An)/dT vs. T and the q = peak intensity vs. T for the d = 2 RFIM system Rb2Coo.s5Mgo.15Fi. 
The birefringence data show that the application of the random field destroys the transition. The data show 
no hysteresis near T C (H), which indicates equilibrium behavior. The neutron scattering peak intensities are 
obtained after ZFC and FC. At low temperatures the long-range antiferromagnetic order Bragg component 
is stable. As T is increase, the long-range order becomes unstable and decays, well below the temperature 
region of the destroyed phase transition. No long-range order is observed upon FC. 



in Rb2Co x Mgi- x F^ by the birefringence experiments of Ferreira, et al. £3 and (dM/dT)n 
experiments of IkedaEl. -^55^ is proportional to the magnetic specific heat^l and is particu- 
larly important for low dimensional systems where the phonon specific heat is considerable. 
The transition for H = is well described experimentally by a symmetric logarithmic di- 
vergence. However, even relatively small applied fields round the transition, as is evident 
in Fig. 5. This behavior contrasts greatly with d = 3 rounding observed upon FC since 
the d = 2 crystal is in equilibrium above and below the H = transition and no hystere- 
sis is observed upon FC and ZFC. As the field increases, the rounding also increases in a 
way consistent with the random-field scaling function (Eq. |i~2|) with a crossover exponent 
<f> = 1.75, which is approximately equal to the zero-field staggered susceptibility exponentS 
as expected. 

Whereas the behavior near TV is in excellent accord with theory, the physics of d = 2 di- 
lute antiferromagnets at low temperatures may not be equivalent to that of the ferromagnet 
with random fields 00. The low T behavior is discussed in the next section. 

The neutron scattering line shapes were studied when the sample was FC to temper- 
atures well below the destroyed phase transition EH where nonequilibrium behavior domi- 
nates. The Lorentzian plus squared-Lorentzian line shape of Eq. |l7] works quite well as does 
a Lorentzian to a power of approximately 3/2. Although no comprehensive study has been 
made of the line shapes near T C (H), preliminary indications are that the squared-Lorentzian 
scattering term is relatively unimportant in Rb2Coo. 85^50.15-^4 111 this equilibrium regionB. 
This conflicts with the mean-field theory that predicts that the squared-Lorentzian should 
be just as important for d = 2 and d = 3 near T C (H) and suggests that the mean-field 
arguments for the squared-Lorentzian are not particularly relevant. This problem deserves 
further study. 
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7 Low Temperature Dynamics in d = 3 and d = 2 

The dynamics of the RFIM below T C {H) in dilute antiferromagnets have been explored using 
Squid magnetometry in Feo.4eZno.54F2 by Lederman, et al.tS After inducing domains using 
the FC procedure, the relaxation of the metastable domain walls has been measured as a 
function of time. The excess magnetization from the domain walls scales as the inverse 
of the domain sizeEi The dynamics for a variety of fields and temperatures have been 
characterized. Below T2(H), which is approximately equal to the equilibrium line T eq (H), 
and above the another line Ti(if), the time dependence of the domain wall size R(t) is 
consistent with the expression introduced by Villain R(t) ~ H~ UH ln(i/r), where r is 
a spin-flip time. This indicates that the dynamics are governed by the pinning from the 
random-field fluctuations. Below T\(H) at lower fields, the random-field pinning seems to 
be insignificant relative to the p inning from vacancies, which are known to freeze in domain 
structure even at zero field Eloi for d = 3. At very low T for all fields the Ising character 
of the spins is sufficient to freeze the domain structure. It is not yet clear how this picture 
might change with variation in the magnetic concentration. The time dependence observed 
by Lederman, et al. is consistent with recent domain growth near T C {H) observed by Feng, 
et al. e3 in Teo.5Zno.5T2 in a very largely field H = 5.5T. The line shape width decreases 
with time near T C {H) but not at low temperatures. The smaller field behavior h as n ot been 
probed. RFIM dynamics have also been observed using Monte Carlo techniques 

The low T dynamics of the d = 2 RFIM dilute antiferromagnet Rb2Coo.s5Mgo.15F4 were 
probed using neutron scattering techniques £3. No Bragg neak develops upon cooling with 
H > since the equilibrium phase transition is destroyedtj. Instead, a non-Lorentzian-like 
scattering line shape develops £3 well below the rounded transition. On the other hand, if 
the system is ZFC, long-range order is observed to be stable at low T for H > 0. Upon 
heating, a temperature region is reached where the Bragg scattering peak decays. This 
region, shown in Fig. 5, is well below the destroyed phase transition as seen by comparing 
with the birefringence data @, also in Fig. 5. The time dependence of the decay of the 
Bragg intensity at the steepest slope in the Bragg intensity, Tp, versus T is observed to be 
approximately logarithmic. Furthermore, the scaling behavior T/v — Tp ~ H 2 ^ is observed 
with cj) = 1.74 ± 0.02, in good agreement with the random-field crossover exponentH3 <fi rj 
1.75. Hence, the instability of the long-range order is certainly connected with the random- 
field behavior. Just as in the case of d = 3, once the domains are introduced into the 
system below T/v and the field is turned off, the domains remain for T < Tv even though 
the ground state is long-range order. 

The dynamics of domain formation at low temperatures have been studied very close 
to the percolation threshold in Rb2Coo.mMgo.40F4 by Ikeda, et al. Il^ using neutron and 
magnetization techniques. Currently, the behavior is being investigated!^ at higher con- 
centration in Rb2Coo.sbMgo.15F4. 



8 Griffiths-like Phase in Dilute Antiferromagnets 



Griffiths LLHil showed that the magnetization in dilute magnets is nonanalytic in H at H = 



below the transition temperature of the corresponding pure system. This is a consequence of 
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the randomness of the local magnetic concentration. Evidence for dilution-induced Griffiths 
instabilities has been observed t3 by studying the deviations from the Curie- Weiss behav- 
ior of x' which appears at the pure Neel temperature and extends down to the transition 
temperature in Feo.47Zno.53F2 and K2Cuo.sZno.2F2. A similar, but much stronger effect 
is observed in Feo.47Zno.53F2 once random fields are introduced. Binek and KleemannEH 
were able to describe the field-induced Griffiths- like peak in seen as the broad peak in 
Fig. 4, using a phenomenological Lorentzian density distribution of local critical tempera- 
tures between T c (0) = T/v and T C {H) with a corresponding power lawv' behavior at each 
temperature. These phenomena have only recently been investigated 113 in d = 2 systems. 



9 The d = 3 RFIM at Large Magnetic Dilution and Large Fields 

New physics emerges once the percolation threshold x ~ 0.24 in Fe x Zn\^ x Fi is approached. 



The system behaves much like a spin-glass I u nr ua l, as was first discovered by Montenegro, 



et al. P^H 107 ! This behavior takes place even though the frustrating exchange interactions 
in Fe x Zn\- x F2 are very smaller. Near the percolation threshold, even tiny frustrating 
interactions are predicted to become important 0. For Ising systems, it is also expected 
that the dynamics even in zero field should be extremely slow 1122. Both of these may 
contribute to the spin-glass-like behavior, although computer simulations seem to indicate 
that the small frustrating interactions are sufficient OU. Very close to the percolation 
threshold, for x = 0.25 and x = 0.27, no Bragg peak, arid hence no antiferromagnetic 



ordering, is observed in zero field with neutron scattering 1123. (Interestingly, this does not 
seem to have been observedjri the related anisotropic systemll2!a Coo.2eZno.74F2 or in the 
weakly anisotropic system 1113 Mn x Zn\- x F2-) The antiferromagnetic correlation length £ 
remains small and constant for T below approximately 10 K (T/v = 78.4 K for pure -Fei^)- 
Mossba,uer measurements indicate a competition between antiferromagnet and spin-glass- 



like order tiH. The temperature below which £ remains constant is just the endpoint of the 



de Almeida-Thouless line T eq (H). The T C (H) curvature is described well 
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ay a crossover 
For a higher 



exponent 4> = 3.4, the same exponent measured in canonical spin-glasses 
concentration, x = 0.31, a more complicated phase diagram is observed 112^. The low- 
field behavior is the same as observed for higher concentrations, i.e. the low-field phase is 
antiferromagnetic and <j> = 1.42. As the field increases, the curvature changes to 4> = 3.4 
and no antiferromagnetic order is observed below T eq (H). The large field induces the 
spin-glass-like behavior away from percolation. As we move to even higher concentrations, 
x = 0.5, very high fields are needed to probe the region above the antiferromagnetic phase, 
as shown by Lima, et al. ti3 employing high-field magnetization measurements. Computer 
simulations Ha indicate that below x ~ 0.6, weak frustration affects the ordered state of the 
REIM in dilute antiferromagnets. 

In the less an isotropic sys tem Mno. 35 Ztiq. Q5F2, somewhat similar behavior to that in 
Feo, 31 Zno. 69-^2 is observedlH3 in magnetization and ac susceptibility studies. There is some 
indication that the phase diagrams may differ in some respects and this is currently under 
investigation. A de Almeida-Thouless line with cj) = 3.4 is observed for Mn$. 35^^0.65-^2 ■ 

A spin-glass-like phase has also been observed above the mixed phase in Fe x Mg\- x Cl2 
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for relatively large magnetic concentrations^^. Slow dynamics are observed 



r the metastable 
. The memory 



domain structure within the mixed antiferromagnetic-paramagnetic phasell 
of domain structure is preserved upon decreasing the field to zero and even upon field 
reversal. The memory effect is also observed after entering the spin-glass-like phase. 



10 First-order to Second-order Transition in Fe x Mg\^ x Cli 



Recently the metamagneticjtransition in Fe x Mg\- x Cl2 for has been studied optically and 
with computer simulations EJJ. Rounding of the metamagnetic transition is interpreted as 
the driving of the transition from first-order to second-order by random fields and random- 
field-induced domain structure. The domain structure is optically observed to be greatly 
altered by the dilution-induced random fields. This is in accord with predictions that 
quenched impurities Ell and random fields — can drive a phase transition from first-order 
to second-order. The concentration at which the metamagnetic transition becomes second- 
order is estimated to be x = 0.6. For sufficient dilution the first-order nature of the transi- 
tion is lost when the avalanche of domain flipping no longer involve^ infinite length scales. 
Universal behavior is predicted for this nonequilibrium transition 
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11 Other RFIM Systems 



Although a great deal of the experiments shedding light on the RFIM have been done 

has 
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on dilute antiferromagnets, other systems have been studied as well. Kleemann 
reviewed random- field domain states in ferroelectric and structural phase transitions. The 
critical behavior of the RFIM structural phase transition in DyAs x V\- x O^ has been studied 



extensively E±3 and compared to the dilute antiferromagnet. Neutron and light scattering 



experiments have been done on binary mixtures in silica gels Certainly more RFIM 
realizations will be studied in the future and will significantly add to our understanding as 
well as incorporate aspects of the dilute antiferromagnet results. 



12 Conclusions 



There is good reason to be optimistic about achieving a good characterization of d = 3 
RFIM critical behavior in the near future. Experiments are nearly at the point where se- 
rious comparisons between theory and experiment can be made. This is possible since the 
high concentration crystals show no evidence for the formation of domain structure or for 
hysteresis in the line shapes well below the transition, two aspects of the experiments at 
lower concentration that have been severe impediments. It would be interesting to investi- 
gate if the remarkable difference in the behavior at high and low magnetic concentration is a 
result of a concentration critical point below which the long-range order becomes unstable. 
One outstanding problem is the lack of a theoretically derived line shape to use in analyzing 
data; the mean-field arguments are clearly inadequate. When such a theory is developed, 
more reliable critical exponents and amplitude ratios will be derived from the scattering 
data. 
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Two kinds of hysteresis can now be distinguished for d = 3. At low magnetic concen- 
trations, vacancies cause irreversibilities and domain formation which are most evident in 
scattering experiments. At all concentrations where transitions take place, there appears 
to be hysteresis, observable in all experiments, that may be attributable to random-field 
critical dynamics. For d = 2, hysteresis occurs only at low temperatures, well below the 
rounded transition. The dynamics of domain formation in this region are still being studied. 

Griffiths-like domain structure dominates the ac suseptibility in the d = 3 random-field 
transition in Feo. 46-Zno. 54-^2- It remains a task to determine whether the small critical peak, 
recently discovered, yields power-law or activated dynamics. 

Near the percolation threshold, it appears that the d = 3 Ising system Fe x Zn\_ x F2 
behaves very much like a spin-glass despite having only small frustrating interactions. 
The behavior in the more isotropic Mn x Zn\- x F2 is being studied to elucidate the role 
of anisotropy in the spin-glass-like behavior. In related studies, intermediate concentration 
crystals of Fe x Zn\- x F2 are being studied in the high-field limit. 

Recent experiments have addressed the random-field effects on first-order transitions 
in Fe x Mgi_ x Cl2- The first-order transition appears to be driven to be second-order with 
sufficiently strong random fields, in agreement with theory. 

A reasonable understanding of the random-field Ising model as realized in dilute anti- 
ferromagnets is emerging, though there is considerable work yet to be done. All aspects 
of the rich behavior of these dilute antiferromagnets are important to characterize partly 
for their intrinsically interesting properties and partly because other materials may show 
one or more of the characteristics. The antiferromagnets are the best studied and probably 
the most easily understood systems. In trying to understand the behavior in more com- 
plex systems, one will have to keep in mind the array of possible behaviors. Certainly, the 
random-field physics will be incorporated into the descriptions of many important materials 
in the future. 
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